Quasi-Nambu-Goldstone Modes in Bose-Einstein Condensates 
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We show that quasi-Nambu-Goldstone (NG) modes, which play prominent roles in high energy 
physics but have been elusive experimentally, can be realized with atomic Bose-Einstein conden- 
sates. The quasi-NG modes emerge when the symmetry of a ground state is larger than that of the 
Hamiltonian. When they appear, the conventional vacuum manifold should be enlarged. Conse- 
quently topological defects that are stable within the conventional vacuum manifold become unstable 
and decay by emitting the quasi-NG modes. Contrary to conventional wisdom, however, we show 
that the topological defects are stabilized by quantum fluctuations that make the quasi-NG modes 
massive, thereby suppressing their emission. 
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Introduction. — Spontaneous symmetry breaking gen- 
erally yields Nambu-Goldstone (NG) modes, which play 
a crucial role in determining low-energy behaviors of 
various systems from condensed matter to high energy 
physics. Although the NG theorem guarantees that NG 
modes do not acquire mass at any order of quantum cor- 
rections, we sometimes encounter situations in which soft 
modes appear which are massless in the zeroth order (tree 
approximation) but become massive due to quantum cor- 
rections. They are called quasi-NG modes pQ and were 
introduced in the context of gauge theories with sym- 
metry breaking by Weinberg • It was shown that such 
modes emerge if the symmetry of an effective potential of 
the zeroth order is higher than that of gauge symmetry, 
and the idea was invoked to account for the emergence of 
low-mass particles. On the other hand, Georgi and Pais 
demonstrated that quasi-NG modes also occur in cases 
in which the symmetry of the ground state is higher than 
that of the Hamiltonian [3] . 

Later on, quasi-NG modes turned out to play promi- 
nent roles in technicolor g] and supersymmetry [S] , both 
of which were proposed as candidates beyond the stan- 
dard model for unification of fundamental forces. The 
quasi-NG modes of the Weinberg type [2] become an im- 
portant ingredient in the physics of technicolor, which 
is a model to avoid a hierarchy problem because it does 
not assume the Higgs scalar field as an elementary par- 
ticle. In this model, the quasi-NG mode is related to the 
vacuum that is energetically selected from among a large 
degenerate family of vacua (vacuum alignment problem) . 
On the other hand, the quasi-NG modes of the Georgi- 
Pais type [5] are inevitable in theories with supersymme- 
try, which is among the most powerful guiding principles 
in contemporary elementary particle physics. This is be- 
cause the vacuum condition in supersymmetric theories 
invariant under a group G is always invariant under its 
complex extension G c . This type of the quasi-NG mode 



TABLE I: Vacuum manifold M and the enlarged vacuum 
manifold M of the uniaxial nematic (77 = rnr/3), biaxial ne- 
matic [77 = (n + l/2)7r/3], and dihedral-2 (other 77) phases, 
where the parameter 77 characterizes the order parameter of 
the nematic phase [see Eq. pb]. 



Phase 


M ' 


* G/H 
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Uniaxial nematic 


U(l) 


x S 2 /Z 2 






Biaxial nematic 


[U(l) x 


SO(3)]/D 4 


[U(l) 


x S 4 ]/Z 2 


Dihedral-2 


U(l) x 


SO(3)/D 2 







is also believed to appear in the weak-coupling limit of 
A-phases of 3 He [B] and spin-1 color superconductivity 
7\. Despite their importance as described above, the di- 
rect experimental confirmation of the quasi-NG modes 
has yet to be made. 

In this Letter, we point out that a spinor Bose-Einstein 
condensate is an ideal system to study the physics of 
quasi-NG modes. This system has recently been a sub- 
ject of active research because of the great experimental 
manipulability and well-established microscopic Hamilto- 
nians. We show that quasi-NG modes appear in a spin-2 
nematic phase, which may be realized by using 87 Rb and 
a cZ-wave superfluid. In the nematic condensate, three 
phases, each of which has a different symmetry, are en- 
ergetically degenerate to the zeroth order [5] as listed in 
Table [TJ We point out that this corresponds to the vac- 
uum alignment problem in particle physics, and to fully 
understand this, we must consider the degrees of free- 
dom of quasi-NG modes due to the U(l) x SO(5) sym- 
metry to the zeroth order. Considering this symmetry, 
we show that the vacuum (order parameter) manifold is 
enlarged to M ^ [U(l) x S 4 ]/Z 2 . Then, the vacuum 
alignment problem can be understood from the fact that 
each vacuum manifold of the NG modes M = G/H is 
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FIG. 1: (color online). Decay dynamics of a Z2 vortex that is 
stable under M = U(l) x S 2 /Z2 in the uniaxial nematic phase. 
(a)-(d) Time evolution of the absolute square of the spin- 
singlet trio amplitude | A30 1 2 defined in Eq. 1 10 1, (e)-(h) Time 
evolution of the quasi-NG charge Q = {F35 cos 8 + F13 sin 9) , 
where 6 is the polar coordinate on a plane and F13 and F35 
are introduced below Eq. Q. The quantum fluctuations 
that stabilize the vortex state are switched off at t = 0. Then 
the vortex decays in time as shown in (b)-(d). This decay 
is accompanied by the emission of quasi-NG modes as shown 
in (f)-(h). The numerical simulations are performed by using 
a modified Gross-Pitaevskii analysis based on a local density 
approximation |14| . In the simulations, we choose numerical 
parameters for typical experimental situations for spin-2 87 Rb 
condensates with n = 1.5 x 10 14 cm -3 . 



a submanifold of M . We show that the vacuum align- 
ment problem is solved because of the quasi-NG modes 
becoming massive and that one of the vacuum manifolds 
is selected because the symmetry at the zeroth order is 
broken by quantum corrections. We find that a soliton 
that is unstable classically becomes stabilized by quan- 
tum corrections that suppress emissions of the quasi-NG 
modes as shown in Figs, [ija) and (e). Notably, as soon 
as the quantum fluctuations are switched off at t = 0, 
the soliton begins to decay since the original vortex core 
structure disappears as shown in Figs, (b)-(d); this is 
accompanied by emissions of quasi-NG modes as shown 
in Figs, (f)-(h). 

Symmetry of the ground state. — We start with a spin- 
2 condensate with mass M whose G = U(l) x SO(3) 
invariant Hamiltonian is given by 



H = 



ci 



2M 



4 A i 



(V>. 



CO . 

2 



(1) 



where ip m (m = 2, 1, — 2) are field operators with mag- 
netic quantum number m, i* 1 , (i = x, y, z) are spin ma- 
trices, T denotes the time reversal operator [Tip m = 
(— l) m i/j_ m ], and repeated indices are assumed to be 
summed over 2,1,...,— 2. The coupling constants are 
related to the s-wave scattering lengths as in the to- 
tal spin S channel by cq = 47rfi, 2 (4a2 + 3ai)/7M, c\ = 



47rft 2 (a 4 -a 2 )/7M, and c 2 = 4irti 2 (7a -lQa 2 +3a4)/35M . 
By applying the Gross-Pitaevskii theory to the Hamilto- 
nian ([lj under the assumption that the k = compo- 
nents are macroscopically occupied, we can show that the 
ground state of the nematic phase is realized for c 2 < 
and c 2 < 4ci [HI DZ3] ■ The order parameter of the nematic 
phase is given by 



(2) 



0m = (sin 77/ v 2, 0, cos 77, 0, sin 77 



where 77 is an additional parameter independent of G [H] , 
and 77 = 7i7r/3, 77 = (n + 1/2)tt/3, and the other values 
correspond to the uniaxial nematic, biaxial nematic, and 
dihedral-2 phases, respectively, where n 6 Z. Since the 
state with 77 is equivalent to that with 77 + 77/3 up to 
U(l)xSO(3), we restrict the domain of definition to < 
77 < tt/3 unless otherwise stated. 

Next, we show that the zeroth-order solution of the 
nematic phase has G =U(l)xSO(5) symmetry. As dis- 
cussed in Ref. [H] , the interactions proportional to Co, c% , 
and c 2 have SU(5), SO(3), and SO(5) symmetries, respec- 
tively, in addition to the U(l) symmetry. On the other 
hand, the zeroth-order solution for the nematic phase 
satisfies |^T^ m | = 1 and cj)l n (F) mn (j) n = 0. In fact, 
the first relation is sufficient to characterize the nematic 
phase, since it is proved that \{ij)] n Tip m )\ = 1 implies 
(F) = [TO]. Within the zeroth order, the interaction 
term proportional to c 2 satisfies |^T0 m | 2 = 1, and is 
invariant under a G transformation, which implies that 
I^J^T^ml = 1 is preserved under the same transforma- 
tion. Therefore, we conclude that the solution of the 
nematic phase has the G symmetry even if c\ ^ 0. 

Number of NG and quasi-NG modes. — We identify the 
enlarged vacuum manifold of the nematic phase based on 
G and discuss the number of NG and quasi-NG modes 
in each phase. Since the vacuum manifold is indepen- 
dent of the location in the orbit, we choose the config- 
uration <f> m = (— 7,0, 0, 0,i) T /^/2, which describes the 
five-dimensional representation of SO(3) in the spheri- 
cal tensor basis and therefore describes the SO(5) in the 
same basis. To determine the isotropy group, let us next 
transform from this representation of SO(5) to the fun- 
damental representation of SO (5). By using a unitary 
matrix 
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1 



( i -i\ 

-i -i 

y/2 

1 0-10 

V-i 00 0-1/ 



(3) 



we can achieve it for the state and generators as 
^ = [/^ m = (l,0,0,0,0) T 



(4) 



and (xkl)ij = {UF k iW)ij = -i{8ki&lj ~ Suhj), respec- 
tively, where F^i is the corresponding SO (5) generator 
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in the spherical tensor basis. In the fundamental repre- 
sentation, elements of the SO(5) group are expressed as 
Rki(t) = exp(-ixkit) with real parameter t. Then, Q is 
invariant under 7? 23 , R24, R25, R34, R35, and -R45, which 
constitute an S0(4) as a subgroup of SO(5). Therefore, 
the enlarged vacuum manifold except for discrete groups 
is U(l)xSO(5)/SO(4)=U(l)xS 4 . To consider the effect 
of discrete symmetry, let us utilize the degrees of freedom 
of U(l)xS 4 . Then, cp/. can be transformed as follows: 

ip' k = e 1 ^ (cos 0\ , sin 9\ cos 62 , sin 9\ sin 6* 2 cos #3 , 
sin 61 sin 62 sin 9 3 cos # 4 , sin 9\ sin 9 2 sin 6> 3 sin 9i) T , (5) 

where < 6\, 62, 63 < tt and < 9 4 < 2tt. To ensure that 
ip' k coincides with cp/., it is necessary and sufficient that 
6\ = (j) = 0, or #i = cj> = 7T, which is isomorphic to Z 2 . 
The full enlarged vacuum manifold of the nematic phase 
is therefore given by 

- ^ U(l) x S 4 ^ U(l) x SO(5) 

Z 2 Z 2 k SO(4) ' 1 j 

Here, x implies that the nontrivial element of Z 2 does 
not commute with some elements of SO (4). 

Next, we discuss the number of NG and quasi-NG 
modes. As analyzed in Ref. [T3], the number of the NG 
modes in the nematic phase is equal to the dimension of 
M, dim(M). Meanwhile, in Ref. 3J, it has been shown 
that the number of quasi-NG modes, n, is given by 

n = dim(M) - dim(M), (7) 

where M is the surface on which the effective potential 
assumes its minimum value to the zeroth order. This 
implies that M is a submanifold of M and n is the di- 
mension of the complementary space of M inside M. 

In the nematic phase, we expect 5 soft modes because 
dim(M) = 5, which is consistent with the result of the 
Bogoliubov theory since it predicts an equal number of 
massless modes |T2j [13] . While three and two of them 
are the NG and quasi-NG modes for the uniaxial nematic 
phase because of dim(M)=3, four and one of them are the 
NG and quasi-NG modes for the dihedral-2 and biaxial 
nematic phases because of dim(M)=4. (See Table |h| ) 
That is, the number of NG and quasi-NG modes changes 
in each phase, with the total number unchanged. 

Quantum symmetry breaking. — Since the interactions 
proportional to Cq and c 2 favor G but the interaction pro- 
portional to ci breaks G, it is expected that the ground- 
state symmetry of the nematic phase at the zeroth order 
is broken by quantum corrections, thereby making quasi- 
NG modes massive. In fact, the ground-state energy per 
particle at the 1-loop level is evaluated as [HI [T3] 

Ae = v(rj) + v, (8) 
= w£ 2 =0 (l + Acos(2n + 27rj/3)) 5 / 2 , (9) 



TABLE II: Broken generators within (2| and < 77 < tt/3, 
the number of ordinary NG modes [U(l), F x , F y , and F z ], and 
that of quasi-NG modes (F13 and -F35) in each phase, where 
F x = -F14 - F 25 + V3F 23 , F y = -F12 + F45 - v^fk, and 
F z — 2Fi 5 + F24- 



Phase 


Broken generator 


N NG 


A'quaai-NG 


Uniaxial nematic 


U(l),F x ,F y , 


Fi3,F 35 


3 


2 


Biaxial nematic 


U(l),F x ,F y , 


Fz,F 35 


4 


1 


Dihedral-2 


U(l),F w ,Fy, 


F Z ,F 35 


4 


1 



where u = 8\/M 3 [n(2ci - c 2 )} 5/2 /15mr 2 /i 3 , X = 
— 2ci/(2ci — c 2 ), v describes an 77- independent contribu- 
tion, and v(rj) is the G symmetry breaking contribution, 
which favors the uniaxial nematic phase for c\ > and 
the biaxial nematic phase for c\ < 0. In other words, 
the degeneracy at the zeroth order is lifted and one of 
the phases is selected, depending on the sign of C\. This 
phenomenon corresponds to the vacuum alignment. Ad- 
ditionally, since v[rj) breaks the extra flat directions, 
masses of the quasi-NG modes that are of the order of 
lj arise. In this way, the quantum symmetry breaking of 
the U(l) x SO(5) symmetry indeed occurs. 

Fate of topological defects. — Here, let us discuss the 
importance of the quasi-NG modes in terms of topolog- 
ical defects. If the symmetry breaking contribution ^ 
is negligible, some topological defects predicted from M 
become unstable. This is because, under such a situa- 
tion, topological defects are characterized by ir n {M) in- 
stead of 7r„(M), where 7r„ represents the nth homotopy 
group. To illustrate this in detail, by using a modified 
Gross-Pitaevskii equation based on a local density ap- 
proximation [14] . we examine the (in)stability of the Z 2 
vortex in the uniaxial nematic phase, which is predicted 
from 7Ti(S 2 /Z 2 ) = Z 2 . In the Z 2 vortex, the spin-singlet 
trio amplitude 

A 30 = 3^/6(^2^.1 + #0-a)/2 + <M0o ~ 3< ^-i - 60 2 </)_ 2 ) 
= cos3?7 (10) 

becomes zero at the vortex core and has a finite value 
in the other region. Initially, the Z 2 vortex is prepared 
under the condition that the contribution ^ cannot be 
neglected, and, therefore, the vortex remains stable and 
the quasi-NG modes are not emitted spontaneously (Figs. 
[TJa) and (e)). However, when the contribution ^ is neg- 
ligible, the quasi-NG charge expands in a radial direction 
and the vortex core structure disappears with time, which 
indicates that the vortex decays due to the emission of 
the quasi-NG modes (Figs. [TJb)-(d) and (f)-(h)). We 
emphasize that this drastic change in the vortex state is 
caused by quantum fluctuations. 

Geometrical interpretation. — The mathematical struc- 
ture of the nematic phase can be interpreted from the 
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FIG. 2: (color online), (a) Graphical interpretation of the 
NG-quasi-NG change within < r\ < n/3. Note that this is a 
vertical slice about the directions along which the NG modes 
on U(l), F x , and F y fluctuate, (b) A simple example of the 
NG-quasi-NG change in which the target space is the rugby 
ball (homeomorphic to S 2 ), the base space S 2 /U(l), and the 
structure group U(l). Since U(l) symmetry is broken at the 
generic points, we have one NG mode along the U(l) orbit and 
one quasi-NG mode along the base space S 2 /U(l) orthogonal 
to the U(l) orbit. However, two quasi-NG modes appear at 
the singularities of the base space where the U(l) symmetry 
is restored. This is because neither of the two directions is 
generated by U(l) at the U(l) fixed point. 

viewpoint of differential geometry [15] . The target space 
M is parametrized by both NG and quasi-NG modes. 
To identify the quasi-NG modes, we eliminate the NG 
modes by taking a quotient of M by G. Then, the quasi- 
NG mode corresponds to the orbit space M/G and it is 
parametrized by quantities made of vb m , which are G- 
invariant but not G-invariant. There should exist n such 
invariants with n defined in ^ at generic points, where 
the dimension of unbroken symmetry dim(_ff) is the low- 
est. For the nematic phase, where the generic points are 
the dihedral-2 and biaxial nematic phases with n = 1, 
there should be one polynomial to describe M/G. In 
fact, we can show that this is uniquely determined to 
be the absolute value of the spin-singlet trio amplitude 
|j43o|. From this result, we find that the orbit space M/G 
can be characterized by r\. Physically, the direction of 
ij is that of the quasi-NG mode, which is generated by 
F 35 . The target space M can be regarded as a (singular) 
fiber bundle over the base space M/G, with fiber M and 
structure group G. The fiber shrinks at certain points 
where H is enhanced and which correspond to the sin- 
gular points in M/G. Since the number of NG modes is 
reduced at such points, some of the NG modes need to 
change to quasi-NG modes to keep the total number of 
soft modes unchanged. In fact, the NG mode generated 
by F z changes to the quasi-NG mode generated by F i3 
in the uniaxial nematic phase, which corresponds to a 



singularity in M/G, as shown in Fig. [5j 

In conclusion, we have shown that quasi-NG modes 
can be realized with the spin-2 nematic phase because 
the vacuum manifold at the zeroth order is enlarged to 
M and each M is a submanifold of M . This corresponds 
to the vacuum alignment problem, where a unique vac- 
uum is selected by considering quantum fluctuations. We 
expect that the phenomena illustrated in Fig. [I] open up 
a new possibility in quantum vortices that can be tested 
experimentally by controlling quantum fluctuations. 
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